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We propose a new method for obtaining the four-dimensional effective gravitational theory for 
five-dimensional braneworld models with arbitrary numbers of branes in a low energy regime, based 
on a two-lengthscale expansion. The method is straightforward and computationally efficient, and 
is applicable to both compactified and uncompactified models. Particular emphasis is placed on 
the behavior of the radion modes of the model, while the massive effective fields are automatically 
truncated. Generally, the radion modes are found to form a (N-l)-dimensional nonlinear sigma 
model. We illustrate the method for a N-brane model in an uncompactified extra dimension. 



1. INTRODUCTION AND SUMMARY 



The past ten years has seen a flurry of activity related to extra-dimensional models. Motivated largely by string 
theory and M-theory, extra-dimensional models have displayed the potential to provide natural solutions to issues 
such as the hierarchy and cosmological constant problems, as well as providing interesting models for dark matter. 
Braneworld models, motivated by the works of Arkani-Hamed et al. pQ and Randall and Sundrum 2, 3j, have become 
a very active field of research, with many papers investigating extensions to the basic ideas (see, eg, [1HS] and citations 
therein) . 

Several different approximation and computational methods have been used to extract physical predictions from 
extra-dimensional models. In particular, many models have an effective four-dimensional regime at low energies, 
where the radius of curvature of spacetime measured by four-dimensional observers is much larger than a certain 
microphysical lengthscale. We review some of the computational methods that have been used to obtain a four- 
dimcnsional description of five-dimensional braneworld models, in order to place our results in context. 

One method is to linearize the higher dimensional equations of motion about simple background solutions, then 
specialize to the long-lengthscale limit in order to obtain the linearized four-dimensional effective theory (which 
roughly corresponds to discarding the Kaluza-Klein modes). This method was used by Garriga and Tanaka [7] in 
their analysis of the RS-I model [5] , who showed that linearized Einstein gravity is recovered on one of the branes in a 
particular regime. Further analyses to quadratic order and analyses on other backgrounds have also been performed; 
see, for example, Refs. [8tillj. Linearized analyses have many advantages: they are quick and simple, and serve to 
identify all of the dynamical degrees of freedom in the theory, particularly the Kaluza-Klein modes. However, the 
linearized method is inherently limited and cannot describe strong field phenomena such as cosmology and black holes. 

A second method is to project the five-dimensional equations of motion onto a brane; see, for example, Ref. [12] , 
This "covariant curvature" formalism fully incorporates the nonlinearities of the theory. However, the projected 
description includes nonlocal terms, and the truncation to a low energy effective theory is nontrivial, except in cases 
with high degrees of symmetry. 

In order to overcome some of these shortcomings, Kanno and Soda 13 15 suggested a perturbation expansion of 
the covariant curvature formalism known as the "gradient expansion method" , which involves expanding the theory 
in powers of the ratio between a microphysical scale and the four-dimensional curvature lengthscale. This approach 
allows a low-energy description of the model to be found, while retaining the nonlinearities of the theory. This method 
has been particularly successful in investigating the cosmology of braneworld models [131 115H17) , and has the benefit of 
providing an explicit calculation of the five-dimensional metric, but is algebraically complex and requires assumptions 
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on the form of the metric. 

An alternative approach to obtaining a four-dimensional effective action, discussed by Wiseman [18] , focusses on the 
radion mode of the RS-I model. Treating the radion mode as a deflection of the branes, the approach uses a derivative 
expansion to calculate its nonlinear behavior. Although this method nicely captures the nonlincarities of the theory, 
it is highly nontrivial, and "guesses" the four-dimensional effective action, based on the first order equations of motion 
the method finds. 

A final method involves making an ansatz for the form of the five-dimensional metric in terms of four-dimensional 
fields, and integrating over the fifth dimension to obtain a four-dimensional action. Examples of this method in the 
literature include Refs. [T7 ] H9H22] . The benefits of this method are the automatic truncation of the massive Kaluza- 
Klein modes, and the computational efficiency in dealing strictly at the level of the action. The main drawback is 
that the five-dimensional metric ansatz must be found (or guessed) using another method. 

Four-dimensional effective descriptions typically contain moduli fields (radion modes) which describe the distances 
between branes. Often, such modes appear as massless scalar fields which couple to gravity in a Brans-Dicke like 
manner (see, eg, Ref. [IH]). This occurs in the RS-I model of two branes in a compactified bulk with orbifold 
symmetry, for example. In this model, the radion mode must be stabilized by some mechanism (for example, by using 
a bulk scalar field as in the Goldberger-Wise mechanism [53]), or else the theory is ruled out for observers on the TeV 
brane (see, eg, Refs. [7] [21]). In theories including multiple branes, one expects several radion modes which may have 
nontrivial couplings to one another and to the four-dimensional metric at the nonlinear level. 

One common extension of the RS models is to consider models with more than one or two branes. A variety of 
papers have considered three-brane models, usually on an orbifold (see, eg, [TBI, |2"TI |2"5H2"7] ) . Some special cases have 
been considered for arbitrary TV- brane models, mostly to investigate their cosmological properties [551 I2E]- A few 
papers comment that their methods should extend to arbitrary A-brane situations (eg, [29.), but little analysis has 
actually been performed in this regard. 

In this paper, we present a new method to obtain a four-dimensional effective theory from an A-brane model in 
five dimensions. We assume that matter is confined to branes with the only bulk field being gravity, and we do not 
invoke mechanisms to stabilize the radion modes. The method utilizes a two-lengthscale expansion to find solutions 
to the five-dimensional equations of motion in a low energy regime. We do not require assumptions about the form 
of the metric, or the existence of Gaussian normal coordinates. The method is computationally efficient, and does 
not require the explicit use of the five-dimensional Einstein equations or the Israel junction conditions. Instead, one 
always works at the level of the action. Furthermore, our method is very general and can be applied to various 
models. The method has similarities to the gradient expansion method (see especially [21]), but is computationally 
much simpler, and can deal with multiple branes in a straightforward manner. A particular strength of the method is 
that it performs a rigorous treatment of all radion modes, and automatically truncates massive modes. We present a 
brief example of the method for the case of the RS-I model [2] , before illustrating the method in detail for the case of 
A four-dimensional branes in an uncompactified extra dimension, deriving the four-dimensional effective action for a 
general configuration. 



A. Applicable Models 



In this section, we introduce the parameters, metrics, and coordinate systems used to define the braneworld models 
which we apply our method to. The most basic model assumes that the extra dimension is infinite and not compact- 
ified, but the generalization to circularly compactified and orbifolded systems is straightforward, and is described in 
Section [2 

We consider a system of N four-dimensional branes in a five-dimensional universe with one temporal dimension, 
with coordinates x r = (x°, . . . , a; 4 ). We denote the bulk metric by g rs (a; ), and the associated five-dimensional Ricci 
scalar by i?*- 5 -* . For simplicity, we assume that there are no physical singularities in the spacetime. 

The N branes are labeled by an index n — 0, 1, . . . , N — 1, so that adjacent branes are labeled by successive values 
of n. We assume that the branes are non-intersecting. Denote the n th brane by B n . On B n , we introduce a coordinate 
system w„ = (u>^, . . . , w^). The location of the branes in the five-dimensional spacetime is described by N embedding 
functions n x r (w^). From these embedding functions, we can calculate the induced metric n h ab on B n , 



Xb«) = 



d' 



-5rs 



(1.1) 



We associate a non-zero brane tension a n with each brane B n , and we also take there to be matter fields n 4>(w^) 
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FIG. 1: An illustration of the model a) before and b) after gauge fixing. The bulk cosmological constants, brane tensions, and 
metrics are labeled. 

which live on B n , with their own matter action n S m [ n h ab , n <j)]. 

In between each brane there exists a bulk region of spacetime, which we will denote TZq, ...,1Zn, with lZ n lying 
between branes n — 1 and n. The first (last) bulk region will describe the region between the first (last) brane and 
spatial infinity in the bulk. In each bulk region lZ n we allow for a bulk cosmological constant A n (see Ref. 28 for a 
possible microphysical origin for such piecewise constant cosmological constants). 

Finally, the action for the model is 

t /d(5) \ w-i . N-l 

s[ 5rs ,"* r ,>] = /dV=ff (or <r-A(* r ) )-£*»/ rfW^+^XIX/^L (1.2) 

J V 5 / „ =0 n=0 

where k§ is the five-dimensional Newton's constant, and A(x T ) takes the value A n in TZ n . 



B. Overview of the Method and Results 

Our method works in five steps. 

Step 1: Gauge specialize. From the general action (Eq. ( |1.2[ ) in the model we discuss here), we perform a gauge 
transformation to specialize the metric to the straight gauge |10j . illustrated in Figure [l] 

Step 2: Separate lengthscales in the action. There are two characteristic lengthscales in the model. The first, 
which we call the microphysical lengthscale, is the lengthscale associated with the bulk cosmological constants, which 
is typically assumed to be on the order of the micron scale or smaller. The second lengthscale is the four-dimensional 
radius of curvature felt on the branes. When the ratio of the microphysical lengthscale to the four-dimensional radius 
of curvature is small (the low energy limit), the dynamics of the extra dimension effectively decouples from the four- 
dimensional dynamics, leading to a four-dimensional effective theory. We introduce a small parameter to tune this 
ratio, and use this parameter to perform a two lengthscale expansion of the action. 

Step 3: Solve equations of motion. The equations of motion at zeroth order in this small parameter are 
calculated and explicitly solved. As expected in this type of model, all of the bulk cosmological constants must be 
negative, and at this order, the brane tensions are required to be tuned to a specific value 1 [2 in order to avoid 
an effective cosmological constant on the branes. The solution to the zeroth order equations of motion provides a 
background metric solution, which is perturbed at the next order in our small parameter (the metric is an exact 
solution if the four-dimensional space is flat). 

Step 4: Integrate five-dimensional dynamics. The five-dimensional dynamics of the theory are integrated 
out by substituting the metric ansatz into the action, and integrating over the extra dimension. The action to zeroth 
order in the small parameter is minimized by the solution, leaving only the four-dimensional terms in the action. 

Step 5: Redefine fields. The final step is to redefine fields to cast the four-dimensional effective action in the 
form of a four-dimensional multiscalar-tensor theory in a nonlinear sigma model. In the Einstein conformal frame, 



We consider small deviations from this value in Section 
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the general form of the four-dimensional effective action is given by 

S[g ab , = / d 4 x^g~ -^R[g ab ] - -j AB (<S> C )V a <Z> A \7 b <S> B g ab + ]T n S m [e 2a ^ C )g abl n cj> (f .3) 

^ K 4 J n=0 

where f < A < N — f , arc massless scalar fields (radion modes) which encode the interbrane distances. Also, 
n\ (= 8itGn) is the effective four-dimensional Newton's constant, which is a function of n\ and the bulk cosmological 
constants. Finally, Jab^ ) is the field space metric of the nonlinear sigma model, and a n (& c ) are the brane coupling 
functions. The functional form of both of these depends on the specifics of the model. 

One of the features of the method used here is that five-dimensional gravitational perturbations which give rise to 
massive four-dimensional fields are automatically truncated. The mass scales of these fields are typically of order h/C, 
where C is the microphysical lengthscale of the theory. However, D amour and Kogan [551 123 US] have shown that it 
is possible to have graviton Kaluza-Klein modes where masses are of order h/Cexp(—l/C), where I is an inter-brane 
separation. Because of the exponential factor, these second graviton modes can be ultralight and observationally 
relevant. Although the models we consider are likely to contain such ultralight graviton modes, our method excludes 
their possible contributions to a four-dimensional effective theory. 

Our method has similarities to the gradient expansion method of Kanno and Soda [13H15] . Our small expansion 
parameter coincides with theirs, and the zeroth order solutions from both methods agree in cases where both methods 
are applicable. However, beyond this point, the methods diverge. Our method Taylor expands the action, but not 
the metric as in the covariant curvature formalism. Although higher order corrections to the metric do exist, they are 
intrinsically five-dimensional interactions which are unnecessary for the construction of a four-dimensional effective 
theory, and their contributions to the effective theory are exponentially suppressed within the low energy regime. 
Furthermore, our method arrives at a four-dimensional effective action, rather than working only at the level of the 
equations of motion. This provides for computational efficiency and a more intuitive understanding of the final result. 

This paper is structured as follows. We begin by briefly illustrating our method in the context of the RS-I model 
in Section [2j We then move on to detailing our method in the context of an uncompactified ./V-brane model with an 
infinite extra dimension. In Section [3j we specialize our coordinate systems to a convenient gauge, and in Section [4] 
we discuss our two-lengthscale expansion in detail. In sections [5j and [6j we solve the high energy dynamics to obtain 
the background metric ansatz, and derive the four-dimensional effective action of the model. We then specialize to 
one-brane and two-brane cases, before deriving the general iV-brane action (Section RPfc. Finally, we discuss the nature 
of the approximations employed, before summarizing and concluding in Section [HJ 



2. APPLICATION OF THE METHOD TO THE RANDALL SUNDRUM MODEL 



To briefly illustrate an application of our method, we apply it to the well known case of the Randall Sundrum (RS-I) 
model with a general background. The derivation of results in this section follows the details on the uncompactified 
model treated in the remainder of this paper closely. 

Many papers have used a metric ansatz for the RS-I model (e.g. [HHH]), guessing at the form of the five-dimensional 
metric, and using this to compute the four-dimensional effective action. Such metrics are typically of the form 
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ds 2 = e^ e ^ 7ab (x c )dx a dx b + ( X,y( 2 8 fc C ' y) ) dy 2 (2.1) 



where k = \f— k|A/6. Rather than guessing at the form of the five-dimensional metric, our method derives a five- 
dimensional metric solution, from which the four-dimensional action is calculated. 

The RS-I model contains two branes on a circular orbifold. We consider the circle of circumference 2L, with the 
branes at y — and y = L, with — L < y < L. We will let the y — brane be the Planck brane, and the y — L brane 
be the TeV brane. Points y and — y are identified. To write this in the language of regions used above, we treat the 
regions — L < y < and < y < L as two distinct regions, but identify fields by 4>(—y) = 4>(y)- 



We now follow the computational steps outlined in Section I B 



Step 1. Write the action in the straight gauge [10]. In this gauge, the general metric is given by 



ds 2 = e xixC ' y) ~f ab (x c , y)dx a dx b + <5> 2 (x c , y)dy 2 



(2.2) 
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where det7 = — 1, and we take $ to be positive. The action is given by 



S = 



d 4 x 



dy ■ 



o+ 



dy 



-L+ 



-1-a)-<t [ d 4 xJ^°h-a L [ d A x\[^h 
H ) Jbo Jb l V 



i?( 5 
2 



+ 1 f d ^ x ^h (o^+ + o^- } + j_ f d ^Jzr h { l k+ + + o Sm [0Kbj + L Sm [Lflab 

K 5 JB K 5 JB L V 



(2.3) 



The indices and L refer to the Planck and TeV branes, respectively. h a f, is the four-dimensional induced metric on 
a brane, and a is the brane tension. K + and K~ are the extrinsic curvature tensors on either side of the branes, and 
S m is the matter action on each brane. 

Step 2. Now, expand the action (2.3| to lowest order in the two-lengthscale expansion detailed in Section |4j The 
action to lowest order in this model is given by 



S 



d 4 x 



dy 



o+ 



dy 



\l ab l b c,yl Cd l d a.,y ~ 5 ( X ,yf ~ ±X,yy + ^§X,y ~ 2«§<I> 2 A 



+ 



dy+ I dy] X(x a ,y) (^=7 - l) + / d 4 xy/^fe 2x ^ 



L+ 

















V * 








- O"0 



2 X (i) 



X^ 



!/=-£+ 



0"L 



(2.4) 



Here, x(0) denotes x(x a ,0), and similarly for x(L). The second line in this action includes a Lagrange multiplier (A) 
to enforce the condition det7 = — 1. 

Step 3. Varying the action (2.4 1 with respect to the three fields \, 7 an d 3>, the following equations of motion are 
obtained. 



_ be 

i 



2y 



The following boundary conditions at the branes are also obtained. 

lab,y{y = 0,L) = 

= 0+) = - ^0o$ 



2 K H 2 A 



(2.5) 
(2.6) 
(2.7) 

(2.8) 
(2.9) 

(2.10) 



We now solve the equations of motion. The solution to (2.6) is given by (in matrix notation) 



j(x a , y) = A(x a ) exp ^B(a; a ) jT $(x a , y') e - 2x ^ a ' y ">dy' 



(2.11) 



where A and B are arbitrary 4x4 real matrix functions of x a , subject to the constraint that 7 is a metric. This can 
be combined with (2.8) to yield B = 0, and so 7 is a function of x a only. The only remaining equation of motion is 
then X % = -2k!$ 2 A/3. Defining k = y/- K lA/6, this gives x,y — ±2fc$. Choos e the nega tive s olution, so that the 
brane at y = corresponds to the Planck brane. The other boundary conditions (2.9) and (2.10) yield 



6k 



6k 

T2 



and a 1, 

which are the well known brane tuning conditions. Combining these solutions, the metric solution is then 

2 

dy 2 . 



ds 2 = e x^,y) lab{x c )dx a dx b + (- X ' y{ f k ' y) J 



(2.12) 



(2.13) 



G 



Step 4. We now have the zeroth-order metric solution, which has solved the five-dimensional dynamics. The next 
step is to use this metric in the original action and integrate over the fifth dimension (c.f. |19J). The zeroth order part 
of the action integrates to exactly zero, while the remainder of the action (the original second order terms) yields the 
following four-dimensional effective action. 



S = 



rf 4 -, 



(l - e xW) 



+ X [ 7 „6, >] + L S m \e^ L hab, L A (2.14) 



The constraint det7 = —1 has been relaxed, instead choosing x(0) = 0. 

Step 5. Transforming to the Einstein frame, let g a b — (l — e*^) y a b, and define exp(x(x a , L) / 2) = 

tanh (^K4if(x a )/\/6). Let k\ — kn\ be the four-dimensional gravitational scale. The action in the Einstein frame 
is then given by 



d^x^J—g 



2k\ 2 



(W)(V a ^) 



Jo 
o „ 



' I wJ 9ab '' °^ 



J,, 



,inh 2 ( W 



(2.15) 



This action corresponds to the four-dimensional effective action arrived at by other means, such as the covariant 
curvature formalism [TH HH • 



For the rest of this paper, we confine our discussions to uncompactified iV-brane models. 



3. THE FIVE-DIMENSIONAL ACTION IN A CONVENIENT GAUGE 



In this section, we take the action (1.2 1 and make coordinate choices to simplify the expression. We also separate 
out contributions due to discontinuities in the connection across branes. We specialize the coordinate system to 
that of the straight gauge [JO], and give the action corresponding to Eq. (1.2) in this gauge. Again, while the 



details presented here are specific to an uncompactified extra dimension, they generalize straightforwardly to the 
other situations described previously. 

In general, the five-dimensional Ricci scalar can have distributional components at the branes, as the metric will 
have a discontinuous first derivative due to the brane tensions. It is convenient to separate these distributional 
components from the continuous parts. It is further convenient to use separate bulk coordinates x v n in each bulk 
region lZ n , rather than using a single global coordinate system. We will therefore have a bulk metric in each region 
lZ n , rather than one global metric. We note that the n th brane will then have two embedding functions: n x^(w^) in 
the coordinates x 1 ^ of lZ n , and n Xn+i( w n) in the coordinates x^ +1 of 1Z n +i. 



Combining these modifications, we can write Eq. (1.2) as 



N r /np(S) \ N ~ 1 i r 



71=0 

N-l 



N-l . N-l 



(3.1) 



where n K + is the trace of the extrinsic curvature tensor of the n th brane in the bulk region lZ n +i, and n K is the 
trace of the extrinsic curvature tensor of the n th brane in the bulk region lZ n , where the normals are always defined 



to be pointing away from the bulk region and towards the brane [see Eqs. (3.14) and (3.15) below]. These terms are 
just the usual Gibbons-Hawking terms [3T)] . 

A note on our conventions. Many functions, coordinates and parameters will be indexed by some index n in this 
paper. For coordinates and parameters, the index will always be in the lower right, e.g, x n . For functions, the 
index will be in the upper left, e.g, n g a p- We use capital Greek letters (T, E, 0) to index five-dimensional tensors in 
arbitrary coordinate systems. When we specialize our coordinate system, we will use lowercase Greek letters (a, /?, 7) 
to index five-dimensional tensors. We use Roman letters (a, b, c) for four-dimensional tensors. The metric g refers to 
a five-dimensional metric, while the metric h refers to a four-dimensional metric. Finally, we use the ( — h + + +) 
metric signature. 
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A. Specializing the Coordinate System 



We begin by specializing the coordinate systems in each bulk region. Denote the coordinates by x^ = (x%,y n ), 
where a indicates one temporal and three spatial dimensions. Without loss of generality, we can choose the coordinates 
such that the branes bounding the region are located at fixed y n . Next, choose the y n coordinates such that the branes 
are located at y n = n — 1 and y n = n. In other words, in the brane embedding functions n x^(w^), 

n - 1 y n {w a n _ l )=n-l, (3.2) 
n y n {w a n )=n. (3.3) 

In this way, the first brane will be located at y = y x = 0, and the last brane located at y N _ j = y N = N — 1. The 
n th bulk region 7Z n then extends from y n = n — 1 to y n = n, with the exceptions of the first and last bulk regions, 
which extend away from the branes to =Foo respectively. 

Next, we use some of the available gauge freedom to remove off-diagonal elements of the metrics. Carena et al. 
[lOj have shown that it is always possible to find a coordinate transformation in lZ n of the form i° — > 6° {x^, y n ) to 
make n g ya — while simultaneously maintaining that the branes be located at y n — n — 1 and y n = n. After such a 
transformation, the metric in lZ n can be written as 

n ds 2 = n lab (x c n ,y n )dx a n dx b n + n &(x c n ,y n )dyl (3.4) 

where the sign of n g yy is known from the signature of the metric. We choose the sign of "$ to be positive. 

The brane positions are now hyperplanes located at y n — integer. It is obvious that only coordinate transformations 
for which y — > g(y) (with no x a dependence) can preserve this form for the hyperplanes. With this condition, only 
coordinate transformations for which x a — > f a (x b ) will preserve the form of the metric. Thus, the remaining gauge 
freedom lies in coordinate transformations of the form x a — > f a (x b ) and y — > g(y) such that the positions of the 
branes are preserved. 

For later simplicity, we choose the following parameterization of the four-dimensional metric n "f ab - In each bulk 
region, let 

n J ab «, Vn) = e"*«*') n j ab (x c n , y n ) (3.5) 

such that the determinant of ™7 af) is constrained to be —1. The function exp(™x) is sometimes called the warp factor. 
The metric in lZ n is then 

n ds 2 = e"*«^> n % h {x c n ,y n )dx a n dxl + n ^(x^y n )dyl (3-6) 



B. Embedding Functions, Coordinate Systems on the Branes, and Induced Metrics 



We now specialize the coordinate system ui° on the n th brane B n . We choose the coordinate system on Bo to 
coincide with the first four coordinates of the bulk coordinate system of TZq, evaluated on the brane. Thus, 

°4K) = ( ^K),°y K)) (3.7a) 
= K,0). (3.7b) 

Now, transform the coordinates in the second bulk region by transforming x® such that 

°^K) = K,0). (3.7c) 

Such a transformation only requires a mapping of the form x 1 — > f a {x\), and so the locations of the branes are 
preserved. Next, choose a coordinate system on B\ such that 

1 ^K) = K,1) (3.7d) 

and continue this process until all branes and bulk regions have related coordinate systems. The coordinate systems 
we acquire have the property that for a point V on B n , we have 

n x r n (T) = n x r n+1 (V). (3.7e) 



Note that while the condition (3.7e| implies that the coordinate patches can be joined continuously from one region 



to another in a straightforward manner, they need not form a global coordinate system as they may not join smoothly 
across the branes. 

From the embedding functions in these coordinate systems we can calculate the induced metric on the branes, using 
Eq. As each brane is adjacent to two bulk regions, there will be two induced metrics, one from each bulk region. 

For B n , the induced metric from lZ n is 



n K b (w c n )=e 



lab(w c n ,n) 



while the induced metric from !Z n +i is 



labium n ) 



(3.8) 



(3.9) 



We will restrict attention to configurations where the two induced metrics coincide (as would be enforced by the first 
Israel junction condition [H]). We then have 



= n Kb«) = n Kb«) 



(3.10) 
(3.11) 



Taking the determinant of this expression and using the fact that the determinants of j a b are constrained to be — 1, 
we find that 



Then by equation (3.111, it follows that 



™xK,n) = " +1 x«,n) 



(3.12) 
(3.13) 



The Action 



Now that we have specialized the coordinate systems for every region and brane in our model, we can rewrite our 
action (3.1| in terms of these coordinates. 



We can evaluate the extrinsic curvature tensor terms as follows. Each brane has two normal vectors, one each from 
the two adjacent bulk regions. We define the normal vectors n n at B n to be the inward pointing normals from lZ n+ i 
and 1Z„. Since the branes are at fixed values of the coordinates y n , this gives 



vr«) 



i 



>$«,n) 



On 



as the normal vector from lZ n and 



g 

• , + 1 $(ix;«,n) Vn+1 



(3.14) 



(3.15) 



as the normal vector from lZ n+ i. The vector "n points to the right of bulk region n towards brane n, while n n + 
points to the left of region n + 1 towards brane n. 

For the extrinsic curvature tensors, we have by definition 



fl("s")a(y) 

^ab\ W n) — « - " <- Vfl n, 



,L a 



dwt, 



n + 
n a 



x n+\ =W niVnA 



Evaluating these using the explicit form of the normals, we have 

1 1 



(3.16) 
(3.17) 

(3.18) 
(3.19) 
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X n+l^ab 



To take the trace of the extrinsic curvature tensor, we contract with the inverse induced metric 

n h ab = e - "* ™7 ab = e~" 

We find 

1$~ 



1$ 



(3.20) 

(3.21) 
(3.22) 



In deriving these equations, we used the fact that n ^ ab n j ab = 0, which follows from det( n 7 ah ) = — 1. 



From Eq. (3.6), the determinant of the fivc-dimcnsional metric can be written as 



(3.23) 



We do not substitute J — n j = 1 at this stage; instead we choose to enforce this at the level of the action by a 
Lagrange multiplier (see Appendix [A|. Using Eqs. (3.21), (3.22), and (3.23), the action (3.1) can be written as 



N 

n=Q Jn n 



n R {S) 



2"x(n) A /Z^ 



7J+1, 



"+1$ 



AT-1 „ N-l 

n=0 n=0 



(3.24) 



4. SEPARATION OF LENGTHSCALES 

We now describe the approximation method, based on a two-lengthscale expansion, which we use to obtain a four- 
dimensional description of the system. We begin by defining the appropriate lengthscales, and then detail how the 
theory simplifies in the regime where the ratio of lengthscales is small. 



A. Two Lengthscales 



There are three groups of parameters in our model: the five-dimensional gravitational scale k§, the brane tensions 
{<J n }j and the bulk cosmological constants {A„}. We assume that all parameters in a group are of the same order of 
magnitude, and so will just consider typical parameters a and A. Working with units in which c = 1, the dimensionality 
of these parameters in terms of mass units M and length units L are [k 2 ] = L 2 /M, [a] = M/L 3 , and [A] = M/L 4 . 

We assume that the dimensionless parameter ct 2 k|/A is roughly of order unity; this will be enforced by the brane 
tuning conditions we derive below [see Eq. (5.17)]. Eliminating k 2 , we can then define a lengthscale by 

C = cr/A (4.1) 



and a mass scale by 



M = rr 4 /A 3 



(4.2) 



For a given configuration, we also define a four- dimensional curvature lengthscale C c (y) on each slice of constant 
y, as follows. We take the minimum of the transverse lengthscale over which the induced metric varies, and the 
transverse lengthscale over which the metric coefficient $ varies. In other words, 



10 



where a, &, . . . denotes an ortho-normal basis of the induced metric, R~t~j is the Riemann tensor of the induced metric, 
and dots denote similar terms with more derivatives. 

Thus, for a given configuration, we have two natural lengthscales: the microphysical lengthscale £ — a / K (the same 
for all configurations), and the macrophysical curvature lengthscale £ c (where the c is intended to denote "curvature"). 

B. Separating the Lengthscales 



We now evaluate the action (3.24) in the low energy regime £ c 3> £, in which the theory admits a four-dimensional 
description. We will find that there is a leading order term of order ~ Ai£, and a subleading term of order ~ 
Ai£(£/£ c ) 2 . Our strategy will be to separate the contributions to the action at each order, minimize the leading 
order piece of the action, and then substitute the general solutions obtained from that minimization into the subleading 
piece of the action. The result will be a four-dimensional action which gives the effective description of the system in 
the low energy regime. 



We write the action ( 3.24 1 as a sum S = S g + S m of a gravitational part S g and a matter part S m , where the matter 



part is the last term in Eq. (3.241 and the gravitational part comprises the remaining terms. 

We first discuss the expansion of the gravitational action S g , which is a functional of a bulk metric g a p and brane 
embedding functions n x T . We define a mapping T e which acts on these variables 

T e :(g al3 , n x T )^(g^, n x r 6 ), (4.4) 

where e > is a dimensionlcss parameter, as follows: (i) We specialize to our chosen gauge, (ii) replace the metric 



(3.6) with the rescaled version 

ds 2 = \e^ x °^% b {x c ,y)dx a dx b + &{x c ,y)dy\ (4.5) 

where indices indicating regions have been suppressed, and (iii) leave the embedding functions in our chosen gauge 
unaltered. We may think of e as a parameter which tunes the ratio of the microphysical lengthscale to the macrophys- 
ical lengthscale. As e is decreased, lengthscales on the brane are inflated, and so £ c increases. Thus, as e decreases, 
so does the ratio £/£ c . In particular, we have 

I). "4 

It is important to note that this e scaling does not map solutions to solutions, but just provides a means of keeping 
track of the dependence on the various lengthscales. 

We can construct a one-parameter family of action functionals by using these rescaled metrics in our original action 



(3.24) 



2. 



S g ,e[9 a0 , n x r ] =e 4 S g [g^, n x T e ]. (4.7) 

We can expand this action in powers of e by 

S g ,e [9 a p, n x r ] = S g>0 [g aP ] + e 2 S g2 [g a p] , (4.8) 

where on the right hand side we omit the dependence on the embedding functions since we have used the gauge 
freedom to fix those. The expansion (|4.8| truncates after 2 terms; there are no higher order terms in e. Note that 



there is no O(e) term, as when the action (4.7) is evaluated, terms of 0(e 2 ) arise from contractions in the Ricci scalar 
using g ab (O(l) terms arise from g vv contractions). Terms of order 0(e) would arise from contractions using g ay , but 
as these components of the metric have been gauge- fixed to zero, they are not present. This can be seen explicitly in 
the expansion of the Ricci scalar ( |A.3[ ). As we tune e — > 0, we move further into the low energy regime, and so we 
identify the zeroth order term as the dominant contribution to the action, and the second order term as the subleading 
term. This provides the separation of lengthscales we desire. 



2 The factor of e 4 in Eq. ( |4.7| ) is inserted for convenience, so that Eq. \4.8\ contains terms of O(l) and 0(e 2 ). This is explicitly shown in 
Section [5] 
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Let us now turn to the matter contribution to the action, S m . We expect the matter action to contribute at 0(e 2 ) 
the same order as the subleading gravitational term. To see this, note that the brane tensions scales as a ~ A4 j C? 
where the scales M and C were defined in Eqs. (4.1 1 and (4.2). The matter action will be roughly S m ~ J p d A 
where p is a 4-dimcnsional energy density. The four-dimensional Newton constant k 2 — 8irG is of order n 
by dimensional analysis [c.f. Eq. (6.11 ) below], and so p will be of order 



x, 
C/M 



M 



2 ri rr2' 

1 A—' ,. A—' r . 



(4.9) 



Taking the ratio p/a now gives 



M/CC 2 C L 2 



M/C 3 



(4.10) 



with a rescaled action S m ,e given by (i) 
(ii) rescaling all fields and dimensional constants with dimensions (mass) r (length) s 



Formally, the scaling (4.10) can b e achieved by replacing the matter action S. 
multiplying by e 4 as in Eq. ( |4.7| 

by factors of e~( r+s ). The expansion of the full action is then 



= S Q + e 2 S; 



e 2 [Sg,2 



(4.11) 



It can be seen that given brane tensions tuned to the bulk cosmological constants, a 2 ~ A/rtf, we require that the 
matter density on a brane should be small, so as not to spoil the tuning. This also yields p <C a, which roughly 
corresponds to the separation of lengthscales condition C <C C c . 

We perform this e scaling separately in each bulk region of the model. The contribution to the action from each 
region will separate into zeroth and second order terms. 



C. The Low Energy Regime 



Now that the contributions to each order have been identified, we can minimize the leading order term in the action, 
S . Once general solutions to the equations of motion have been found, we can use these solutions in the second order 
term in the action. Thus, we solve for the high energy (short lengthscale) dynamics first, and use the solution to this 
as a background solution for the low energy (long lengthscale) dynamics. At this point, we may let e — > 1, and rely 
on the ratio (£/£ c ) 2 being sufficiently small to provide the separation of lengthscales. 

The effect of this separation of lengthscales is to enforce a decoupling of the high energy dynamics from the low 
energy dynamics. We will see below that the equation of motion for the high energy dynamics contains y derivatives, 
but no x a derivatives. The theory at this order thus reduces to a set of uncoupled theories, one along each fiber 
x a = const in the bulk. These theories are coupled together at 0(e 2 ), and thus in the regime of interest, the coupling 
is minimal. After solving the high energy dynamics along these fibres, a four-dimensional effective description of the 
system remains. 

The low energy regime, in which the theory admits a four-dimensional description, is the regime 



Cc > C. 

This regime is also frequently characterized in the literature by the condition 



(4.12) 



(4.13) 



where p is the mass density on a brane and it is a brane tension [c.f. Eq. (4.10) above] . One can interpret the condition 



(4.13) as saying that the mass density on the brane must be sufficiently small that the brane-tuning conditions [Eq. 
(5.17) below which enfor ces a 2 ~ A/k 2 ] are not appreciably modified. However, the condition ( 4.13 1 is less general 
than the condition (4.12 ), and although necessary, is actually insufficient. First, as discussed above, (4.13 1 only applies 



to branes, whereas (4.12) applies at each value of y, including away from the branes. Second, even when the density 



on a given brane vanishes, four-dimensional gravitational waves on that brane can give rise to radii of curvature C c 
that are comparable to C. In this case, the separation of lengthscales will not apply and the four-dimensional effective 
theory will not be valid, despite the fact that the condition (4.13) is satisfied. Curvature associated with the metric 
coefficient $ can also yield similar results. 
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Finally, we discuss a subtlety in our definition of the "low energy regime" . As noted in the previous paragraph, C c 
varies with position in the five-dimensional universe. Our separation of lengthscales will break down when the induced 
metric on any slice of constant y has a radius of curvature C c comparable to that of the microphysical lengthscale C; 
it is insufficient to require that C 3> L on each brane. When this happens, the terms of order e 2 will couple strongly 
to the O(l) terms, and our approximate solutions for the five-dimensional metric will no longer be valid. This will 
gcncrically occur at sufficiently large distances from the branes, as exp(™x) typically grows exponentially small away 
from the branes, and C~ 2 cx exp(— n \)R^ ■ Despite this breakdown, the contribution to the action from these regimes 
is exponentially suppressed by the warp factor, and thus provides only a small deviation from the effective theory. It 
is unlikely that the warp factor can grow without bound after encountering this regime while maintaining a globally 
hyperbolic spacetime. 



THE ACTION TO LOWEST ORDER 



In this section, we calculate and explicitly solve the equations of motion to lowest order in the two lengthscale 
expansion. First, however, we write out the complete, rescaled action showing explicitly the dependence on e. In- 
serting the decomposition (A. 3 1 of the Ricci scalar and the rescaled metric (4.5) into the action ( 3.24[ ) [following the 
prescription of Eq. (4.7)], we obtain 



F 2 X / 1 

[ n % b , "<&, n X , H 4>] = E / d " X n V 73 ^ - fr bn %c,y n Y dn j d a,y ~ - 4"X 

+ 4 ^f n X,y ~ 2 4 n ® 2 K ) + "A(x a , y) l) 

N-l 

+ E / ^/ 2 " ri "'fi 

tTr. J Br, 



mi 



71=0 

N 

E e2 

71 = 





(n X,y 




71+1 




4 




Vn~ n 


71+1$ 


y n +i 



3™$V a V a "x 



; *(V a "x)(Va"x) 



2V Q V a "$ - 2(V Q "x)(Va"3>) + 



N-l 
n=0 



t abi 



(5.1) 



where we include the Lagrange multiplier terms (A.2) discussed in Appendix [Aj and where the factor of e 2 in front of 
the matter action comes from the process described in the previous section. This form explains the choice of the e 4 



factor in Eq. (4.7), and shows the decomposition into O(l) and 0(e ) terms, as claimed in Eq. (4.11) 



From the form of Eq. (5.1 ), we see that we can neglect the last two lines in the limit e — > 0. We can obtain a more 
precise characterization of the domain of validity of this low energy approximation by estimating the ratio between 
the terms dropped and the terms retained. As an example, consider the first term on the 4 th line and the first term 
on the first line. Their ratio is (dropping the 'n' labels) 



-i ab lbc.,yl Cd lda,y 



2x 



[*v; 



(5.2) 



where y is the coordinate lengthscale over which 7^ varies. We recognize the first factor as essentially the Ricci scalar 
of the induced metric e XA f a b, which is of order C~ 2 . We recognize the second factor as the square of the physical 
lengthscale in the y direction over which 7 varies, which is always ~ C 2 (see the explicit solution ( 5.20[ ) below). Thus, 
the ratio is (£/£ c ) 2 , confirming the identification of the low energy regime as L <gC C c . 



A. Varying the Action 



In the action (5.1) at zeroth order in e, we have three fields to vary (in N regions): n ^(x c ,y), n x(x c ,y), and 
™7 ob (a; c , y). There is a subtlety in the variation however. The constraint that det ( n 7 ab ) = — 1 must be imposed either 
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at the level of the equations of motion, or by a Lagrange multiplier. The Lagrange multiplier is explicitly included in 
Eq. (5.1 1. Further details are provided in Appendix |A| 

We begin by varying with respect to n <&. From this variation, we find a single equation of motion in each region, 



n 2 



„2 n.i.2 , 



(5.3) 



Next, we vary with respect to n ^f ab - Note that in varying the action, we obtain boundary terms from neighboring 
regions from the relationship (3.13). The variation produces an equation of motion in each bulk region, 



tad, 



<)<! 



iab,y I icd,y 



^ <f> 



(5.4) 



(If using Lagrange multipliers, this equation results after the Lagrange multiplier is eliminated by tracing the equation 
using n j ab and back substituting). Note that tracing over the indices and using Eq. (A.4) leads to Eq. (A. 5) as 
expected. We also find a boundary condition to be satisfied at each brane, 



1 



^lab.yiVn = Tl) 



1 



n+1/ 



«], ■"""-</'-•"' "' ,,+ L<j, lab,y(Vn+l n ) 



(5.5) 



Finally, we vary with respect to n x- Here, we once again have boundary terms arising from integrating bulk terms 
by parts in neighboring regions. There is an equation of motion in each bulk region, 



1 

12 



/ '6c, y J ida,y 



n 3) 

,y n 



0. 



(5.6) 



We also find a boundary condition at each brane, 



1$ 



n+l. 



Vn= n 



(5.7) 



B. Solving the Equations of Motion 



We have three equations of motion for each bulk region, as well as numerous boundary conditions for the fields 
at the branes [Eqs. (3.12|, (3.13), (5.3), (5.4), (5.5), (5.6), and (5.7)]. Note that these equations all describe the 
dynamics along a fibre of constant x a which doesn't couple to any other fibres, and so solving these equations of 
motion consists of solving the dynamics of the extra dimension of the model. 



We begin by solving Eq. (5.4). It is convenient to solve this equation in matrix notation. Let 



[lab 



7 



where we suppress indices n. Then in matrix notation, Eq. (5.4) is 



/■y Ay 



2y 



(5.8) 



(5.9) 



where dots denote derivatives with respect to y. It is straightforward to check that a solution to this differential 
equation in region n is 



7(x Q , y) = A{x a ) exp (B{x a ) £ $(x a , y')e- 2x{xa ' v>) dy' 



(5.10) 



where A and B are arbitrary 4x4 real matrix functions of x a . The lower limit on the integral is chosen so that the 
boundary conditions may be matched at the previous brane (obviously, care must be taken in the first region). The 
expression (5.10) has the correct number of integration constants to satisfy arbitrary boundary conditions. From our 
knowledge of "f a b, A must be a symmetric matrix with determinant — 1. The exponential has unit determinant, and 
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so B must be traceless. The condition that 7 is symmetric implies that B T = A B A 1 . The quantity which appears 
in Eqs. (5.3) and (5.6) is 



I >bc,y I 'da,y ' iab,yy 



= Tr(B 2 (a; a )) $ 2 e 



-4x 



We define 



b(x a ) = -Tr(B 2 (x Q )) 



(5.11) 
(5.12) 



where the factor of 12 has been chosen for later convenience. From combining Eq. (5.5) with Eqs. (3.12) and (3.13), 



we see that B (and thus b(x a )) is independent of region, while A will change with each region according to Eq. (3.13 1 



From Eq. (5.3), we find 



± \/6™$ 2 exp(-4™x) - o K 5™ $2A r, 



P„"<^/6exp(-4«x) ~ 3«|A n 



where P n is either ±1 and is constant in each bulk region. Differentiating Eq. (5.13) gives 



X 



'HI 



(5.13) 



(5.14) 



The same result is obtained by substituting Eq. (5.3) into Eq. (5.6), and so we see that these equations of motion are 
degenerate. This leaves only one equation of motion (Eq. ( |5.13 ijand one boundary condition (Eq. (5.7)) to satisfy. 



Classes of Solutions 



If B(x a ) = 0, then the induced metric on all the branes are related to one another by conformal transformations, 
and a four-dimensional effective action is easily calculated. On the other hand, when B(x a ) 7^ 0, the induced metrics 
on each brane are not simply related conformally, but through the equations (3.13) and (5.10). If solutions with 



B(x a ) ^ were to exist, the four-dimensional effective theory would contain more than one massless tensor degree 
of freedom; i.e., it would constitute a multigravity theory (see Damour and Kogan [55] ). No such degrees of freedom 
have been seen in any linearized analyses 3 . It is important to note that this is not a Kaluza-Klein mode. We believe 
that solutions with B(x a ) ^ are ruled out due to divergences at y — » ±00, leading to a lack of global hyperbolicity in 
the spacetime, although we have been unable to prove this rigorously. We will restrict attention to the case B(x a ) = 
for the remainder of this paper. 



D. General Solutions at Leading Order 



With ~B(x a ) = 0, the field n j ab becomes independent of y [sec Eq. (5.10)], and also independent of n by Eq. (3.13). 
This means that we can drop the index n from x^, w ni and n 7 af) without causing confusion. With b(x a ) = 0, the 
remaining equation of motion and boundary condition simplify somewhat. Equation (5.13) becomes 



X 



,y — V 3' 

which implies that A„ < 0, and so the bulk regions must be slices of anti de-Sitter space. Define 



(5.15) 



(5.16) 



In addition, it can be shown that in orbifoldcd models, there are no solutions with B(x a ) ^ 0; see Section [2] 
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We can use Eq. (5.151 for n \ in Eq. (5.7 1 to obtain 



h P 

^n 1 n 



(5.17) 



These relations are the well-known "brane tunings" , which determine the branes tensions required in order to avoid 
a cosmological constant on the branes [2]. 



We may integrate Eq. (5.15) to find 



'2fc P J? "$(z a ,yW + f(x a ) n = 0, 

^Xix^n- 1) + 2k n P n $l_ x "<f>(x a ,y')dy' n > 



(5.18) 



where f(x a ) is an arbitrary function. Note that the field n x is related to the distance from the previous brane to y 
along a geodesic normal to the branes, made dimensionless by the appropriate lengthscale in the bulk. In particular, 
X describes the number of e-foldings the warp factor in the metric provides between two points in the five-dimensional 
spacetime. Assuming that $ is not divergent, if exp(™x(j/)) approaches zero or oo anywhere, it can only occur as 
y — > ±oo. We will restrict attention to the cases 



P = +1, and P N = -1. 



(5.19) 



When these signs fail to hold, then the warp factor increases monotonically as one goes to infinity, and it seems likely 
that the spacetime cannot be globally hyperbolic. We exclude cases where exp( n x(y)) — > at finite y by restricting 
ourselves to topologically connected spacetimes [2H1 132] ■ 



E. Summary 



We summarize our results so far. We have N branes, each with a brane tension which has been carefully adjusted, 
according to Eq. (5.17|. The branes divide our system into N + 1 regions. Our coordinates are x a , describing 



four-dimensional space, and y, describing the extra dimension. 

We expanded the action in terms of our e scaling parameter to separate the high and low energy contributions. 
Specializing to a low energy regime, we solved for the high energy dynamics, arriving at the metric for each region of 
our system: 



n ds 2 



n " x{x ^ y) %b{x c )dx a dx b 



n X 2 y (x c ,y) 



dy 2 



(5.20) 



with n x given by Eq. (5.181, where ™$(a; a ,y) can be chosen freely. The pa ramet ers k n are determined by the bulk 
cosmological constants and the five-dimensional Newton's constant, by Eq. (5.16). The derivative x.y nas fixed sign 



P n = ±1 in each region, although the derivative may approach zero as y — > ±oo. 



As an aside, when the metric in each region is in the form (5.20) , the zeroth order action So[g a b] [Eq. (4.8)] 
evaluates to exactly zero. This can be seen by substituting the metric (5.20) into the action and explicitly evaluating 



the integral over the y dimension. All of the integrals become total derivatives whose boundary terms exactly cancel 
the boundary terms present in the action at this order. 



The background metric ansatz ( 5.20 ) is essentially the same as the zeroth order metric calculated by Kanno and Soda 
|14j . taking $> 2 (x a ,y) — exp(2cf>(y, x)) in their notation. However, from here, we proceed without their assumption 
that 4>(y,x) — 4>{x). The "naive" ansatz and the CGR ansatz of Chiba [19] are also in the form of our metric (5.20). 



6. THE ACTION TO SECOND ORDER 



In this section, we investigate the action to second order in e. By integrating out the previously determined high 
energy dynamics, we find the four-dimensional effective action. 
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A. Acquiring the Four-Dimensional Effective Action 



Using the metric (5.201 in Eqs. (3.24) and (4.8 1, we can calculate the second order contribution to the action 5*2. 
The result is 



JV 



s 2 [ n % b , n x, n - 



n=0 



E / <^ x nV~^ 



4«; 5 fc n P n 



\ y R {i) ~ 3 n X, v V 2n X ~ \ n X^ an x){V a n x) 



N-l 



2V 2 " X „ " 2(V a " X )(V a "x,,) + E " 5 ™ e " x(x<i ' n) 7a6, n 4> ■ (6.1) 



71=0 



Note that covariant derivatives as written here are associated with the metric jab, as is the four-dimensional Ricci 
scalar RS 4 \ 



To obtain the effective four-dimensional action, we integrate over y in the five-dimensional action (6.1), as the 
dynamics of this dimension have already been solved. The term involving the Ricci scalar can be integrated straight- 
forwardly, as i?( 4 ) has no y dependence, but the other terms require more manipulation. We can combine the last 
four terms in the five-dimensional integral in the following way: 



-3e"* n X ^ 2n X - ^e"* "x„(V Q "x)(Va"x) - 2 e "*V 2 "x,, - 2e"* (V a "x)(V Q "x, y ) 
" 2dy 



(e" x (V an x)(V a " X )) - V Q (3e"* r % y V a n X + 2e"*V a " X , y ) 



(6.2) 



The covariant derivative commutes with the integration over the fifth dimension in the action, and thus gives rise to 
a boundary term at x a — > oo, which we discard. We obtain 



5rn ~ n n 
2 I lab, X, 



N 

E 



i 



d 



n=0 
N-l 



""' r ' v ^APndy 



eXR (4) + |e"*(V°"x)(V "x) 



n=0 



(6.3) 



Integrating over y, we find boundary terms at each brane and at y = ±oo. We note that the integral converges in the 
first and last regions because of the choices Pq = +1 and Pjy = — 1, and so the terms evaluated at ±oo vanish. We 
can rearrange the remaining terms into a sum over the branes. 



N-i r 

S2n ab , n x, n 4>}= £ / d 



1 



x\j —7 



71=0 ' 

N-l 



1 



1 



E " 5 « 



n=0 



4^5 V^n^n ^n+l^n+l 

x(x<, ' n) 7 o6 , V 



^ i? (4) + _ e " X(V a„ x)(Va n x) 



(6.4) 



B. Field Redefinitions 



We now have a four-dimensional Ricci scalar and a number of scalar fields, n x(x a , n), whose values depend on the 
distance from the first brane to the (n + l) th brane. We modify our previous notational conventions as follows. We 
redefine the fields "x and j a b via 



lab 

n x{x a ,y) 



e f{xa) lab{x a ) 

n x(x a ,y)-f(x a ) 



(6.5) 
(6.6) 



where f(x a ) is defined in Eq. (5.18). Thus, we no longer impose the constraint that det(7) = —1; instead we have 
det(7) = — exp(4/). With the relaxation of this constraint, a Lagrange multiplier term in the action to enforce the 
constraint is no longer necessary. As x wm only be evaluated on the branes from now on, we further depart from 
our previous conventions, and write Xn(x a ) to represent n x(x a ,n). Our new conventions also enforce xo — 0, which 
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means that the metric %b is the Jordan frame metric of the first brane, Bq. Note that the five-dimensional metric is 
unchanged by these notational changes. 

For convenience, we define the quantities 



A n 

B n 



1 



V 

1 1 



e„ = sgn 



k n P n fc n ^-i_P n _|_i J 
for < n < N — 1. It is useful to note that e„ can be written as 

e n = - sgn(a„P„P„+i) 



(6.7) 



(6.9) 



(6.10) 



by using Eq. (5.17). We define an effective four-dimensional gravitational constant n\ that is measured by observers 
on the first brane, Bq, by 



1n\ 4k 5 



2 A . 



Finally, we define fields tp n by 



(6.11) 



(6.12) 



for 1 < n < N — 1. The domain of is the positive reals. Rewriting our action in terms of these new variables (and 
suppressing the subscript '2' from now on), we obtain 



S[% b ,iJ) n , n (t>] = / d xy/-*/ 



2k\ 



N-l 



N-l 



i? (4) 1 + £ W£ + 6 £ 



iV-1 



ra=l 



n=l 



"7ab, 



(6.13) 



This is the four-dimensional effective action in the Jordan conformal frame of the first brane, Bq. Note that the target 
space metric, determined by the kinetic energy term for the scalar fields, is flat, and the target space manifold is a 
subset of the quadrant of IR^ -1 in which all the coordinates ifjn are positive, bearing in mind that each i/j n will be 
bounded either above or below by their definition (6.121 and Eq. (5.18). 



C. Transforming to the Einstein Conformal Frame 



The Einstein conformal frame metric for the action (6.13) is g a b = 7 Q &|0|, where is given by 



N-l 



6 = 1+ e 0£nV>n- 



(6.14) 



The four-dimensional effective action becomes 

e sgn(9) 



S[g a b,ip„, n <l 



d x-s/—g 



2k\ 



N-l 



0a 



9abi °<i> 



N-l 
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n=l 



9ab, U (f> 



(6.15) 



where tildes refer to the metric g a b- Note that the kinetic energy terms in this action (6.15) have apparent divergences 
at = 0. However, for any given set of signs e„ (which correspond to a choice of model), it can be shown that |0| is 
bounded away from zero. This arises because of the way each ip n is bounded either above or below. 
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ANALYSIS OF THE ACTION 



In this section, we analyze the four-dimensional effective action (6.15) in a variety of cases. We begin with the cases 



of one and two branes, which serve to highlight some features of the model in the general case. In these special cases, 
our results reduce to previously known results. We then analyze the general situation. 



A. One Brane Case 



In the one brane case, the effective action simplifies greatly. 

or 0,1 1,1 £ „ 

J [9ab, 



A^F)[ s ] + X[ Wl ^]. 



(7.1) 



The four-dimensional effective action is just general relativity (eo = +1 if the brane has positive tension). This 
corresponds to the RS-II model [3]. 



B. Two Brane Case 



Here the parameter of importance is eo£i, which from Eqs. ( 5.17[ ), (5.19) and (6.9) is given by 

e £i = - sgn (er cn) . 



(7.2) 



With Pq and P2 predetermined, it is possible for one brane tension to be negative, but not both. Therefore eo£i is 
positive if there is a negative tension brane, and is negative if both branes have positive tension. Using the definition 
(6.14) of 6, the action (6.15) becomes 



cr nil / ,4 / — e sgn(l + eoe^i) 
S[g a b,ipi, 4>, 4>\ = d x^-g —2 



|1 + e eiV?| 



2l 9ab, V 



R (4) [g}+6 01 (V^XV^r) 
(1 + eoei^l) 

Bijl + eoei^il 



(7.3) 



/. Positive Brane Tensions 



When both branes have positive tension, e e i = — 1- Which of e and t\ is negative depends on the sign of 0. 
Combining Eqs. 16.14b and (|6.12|), 



9 = 1- Sie Xl . 



From Eqs. (5.19) and (6.10), we see that 



-sgn(Pi 



(7.4) 



(7.5) 



Combining this with Eq. (5.18) and recalling that xo = 0, we see that the exponential function in Eq. (7.4) is greater 



than unity for P\ = +1, and less than unity for P\ = —1. If Pi = +1, then the brane tensions [Eq. (5.17)] require 



that fco > and we see that P>i > 1, giving 6 < for e = —1, e\ = +1. If P\ = —1, then the brane tensions dictate 
that fco < fci. Thus, in this case, B\ < 1, and so 6 > for eq = +1, e\ = —1. 
Assuming that < ipi < 1 (6 > 0, Pi = —1, eo = +1), we define 



where 



(f = /itanh 1 {tpi) 



K4 



(7.6) 



(7.7) 
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The domain of <p is to oo. The action (7.3) then becomes 

1 



Oj, 1, 



d X\/^g 
°S 



2K |^ (4) [5]-^(VV)(V a ^) 



cosh 2 | - ) g ab , °- 



1 S 



1 



2 / V 



fix Sinh " U I V 



(7. 



Requiring that the branes do not intersect or overlap gives 



< ipx < y/Bt 



i-h/k 2 
1 + fci/V 



(7.9) 



Note that k\ < k 2 to satisfy Eq. (5.17), and that yBi < 1 (responsible for G > 0). Thus, Eq. (7.9) is a more 
stringent constraint than < ipi < 1. 

In the situation where ipi > 1 (O < 0, P\ = +1, ei = +1), we define 



ip = fx tanh 



(7.10) 



The domain of ip is from to oo. The action (7.3) then becomes 



S[g a b, ipi, °4>, V] = / d 4 x^/^g 
+ °S 



^( 4 )[ 5 ]-^(vv)(v^) 
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sinh 2 ' 



1 S 



1 



2 / f 



(7.11) 



which coincides with the previous action (7.8) if we swap the actions °S m and S m and rescale units in each matter 
action by factors of 



±1/2 



The constraint on the radion field we impose to ensure that the branes do not overlap in this case is 



1 + fcl/fc 2 



1 - fcx/fco 



> 1, 



(7.12) 



where k\ < fco from the brane tunings (Eq. (5.17)) 



The two cases P\ = —1, described by the action (7.8), and Pi = +1, described by the action (7.11), lie on opposite 



sides of 6 = 0. Each set of values eo, . . . , zn-i gives different theories, and each theory comes with its own field space. 

The actions ( |7.8|) and (|7.11 ) coincide with formulae in the literature for the action for the RST model, up to a 
rescaling of units [21 [T9l 120] (also, c.f. Eq. (2.15)). They describe a Brans-Dicke like scalar-tensor theory of gravity, 
with matter on each brane having a different coupling strength to the scalar component. 



77. One Negative Brane Tension 



If e e i = 1 then 6 > always, and by requiring the conditions (5.19), both e and ex must be positive. We define 

<p = /j tan~ 1 (■!/?! ), 



(7.13) 



where the domain of tp is to (tt/2)/x. The action (7.3) becomes 
S[g a b, °<f>, V] = J d 4 x^g 
+ °S, 



2^ 2 



2 I f \ 

cos ( - ] g ab , u - 



l( 3 



1 • 2 ( <P \ 1 

^sm ^)g ab , ■ 



(7.14) 



Note that <p is a ghost field, which gives rise to the usual instability associated with a negative tension brane. 
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C. General Case of N branes 



In the general case of N branes with N > 2, we have N — 1 scalar fields ipi, . . . , ipN-i- It will be useful to regard 
the (N — l)-dimensional space of field configurations as an N — 1-dimensional hypersurface in an iV-dimensional space 
with coordinates tpi, . . . , t/'jv, where 



t/) N = \/-e Q E N Q 



\ 



N-l 



E 6 e n1pr, 



(7.15) 



by Eq. (6. 14 1. We also define cat 



eosgn(6). The action ( |7.3[ ) becomes 
~g e sgn(6) 



^^ (4) [.9] + ^ 2 E^(V^„)(V a ^) 



2k 4 

N-l 



n=l 



2 



BJQ 



gab, n <$> 



(7.16) 



where it is understood that the action is constrained by Eq. (7.15) (this constraint may be enforced by a Lagrange 
multiplier if desired). 

The N-dimensional field space metric (rfS 2 ) is now conformally flat (and indeed, is a form of the metric on N- 
dimensional AdS space with an unusual signature). The physical field space metric is obtained by computing the 
induced (N — l)-dimensional metric on the hypersurface ( 7.15[ ) from the A^-dimensional metric. Of the N — 1 fields 
-01, ■■■ , ipN-i, let P with 0<P<iV — lbe the number of fields with eoe„ positive, and let M — N — 1 — P be the 
number of fields with e e„ negative. We now relabel the fields based on which have positive kinetic coefficient, and 
which have negative kinetic coefficient. Call the fields Pi, 1 < i < P, and raj, 1 < j < M. Following the convention 
in Eq. (1.3) for the metric on the A^-dimensional target space of the scalar fields, we can write it as 



6 



N 

^2 e e n dtpl 

n=l 
P 



A I 



i=l j=l 



sgn{e)dtp 2 N . 



(7.17) 
(7.18) 



As the metric on the positive (negative) coordinates is Euclidean, we can transform them into spherical polar coor- 
dinates. Define new coordinates £, . . . , 0p-\ and r\, \\, . . . , Xm-i, such that 



(pi,-..,pp) = C (cos(0i),sin(0i) cos(0 2 ), 
(mi, . . . , tojv) = rj (cos(Ai), sin(Ai) cos(A 2 ) 



, sin(6>i) sin(6> 2 ) • 
,sin(Ai)sin(A 2 ) 



•sin(^p_i)) , 
• •sin(AM-i)) ■ 



In these coordinates, the field space metric is 

9 



dH z = -d( z - Cdfl z p + chf + ifdVt z m + sgn(9)^, 



(7.19a) 
(7.19b) 



(7.20) 



where dQ 2 — d8 2 + sin 2 (81) d9 2 + ... is the metric on the unit (P— l)-sphere, and similarly for dfl^. All of the angular 
fields (9i and Xj) have a domain of to ir/2, as each tp n is positive. We will refer to this region as the "positive 
quadrant" of the maximally extended fieldspace (where the domains of the angular fields are from to ir or 2ir as per 
usual). The following relationships hold: 



c 2 = Ep? 



i=l 
A I 



e 



1 + C 2 - 



v 2 - 



(7.21) 

(7.22) 
(7.23) 
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The domain of r\ and £ is from to oo. 

We now compute the physical field space metric on the (N — l)-dimensional target space (da 2 ) using the constraint 
(7151). We obtain 



dip 2 N = 



1 



|i + C 2 -^l 

2 



da 



1 + ( 2 - TJ 



(CdC~vdri) 2 , 

l-r, 2 



-dC 



i + C 2 - v 2 



C 2 dn 2 + di] 2 



i + C 2 



1 + C 2 - rf 



ifdtt 2 



2 V ( 



i + C 2 -v 2 



drjdC, 



(7.24) 
(7.25) 



To summarize, our final result for the four dimensional action is a massless multiscalar-tensor theory in a nonlinear 
sigma model, of the form 



1 



1 



B ah 



N-l 



(7.26) 



where we assume that eosgn(0) = +1, which is required for a physical theory, given our conventional choices. The 
scalar fields & A are 

= {C, 0i, . . . , 77, Ai, . . . , A A/ _x} . (7.27) 

The field space metric jab(& C ) is given by <i<7 2 [Eq. (7.25)], and the brane coupling functions a n (Q c ) by 

1 



e 2a„ = 



|i + C 2 -v 2 V 

\l + C 2 - V 2 \B n 



, 1 < n < N - 1, 



where £?„ is given by Eq. (6.8), and ip n is defined by the relevant expression in Eq. (7.19) 



(7.28a) 
(7.28b) 



8. DISCUSSION AND CONCLUSIONS 

We have presented a new method to calculate the four-dimensional effective action for five-dimensional models 
involving N non-intersecting branes in the low energy regime. Although we have only illustrated an application of 
the method to an uncompactified extra dimension, it is generally applicable, and is expected to work for circularly 
compactified and orbifoldcd models also. 



A. Domain of Validity of the Four-Dimensional Description 



We begin our discussion of the domain of validity of the four-dimensional description given by Eq. (6.15) by 
recapping the method of computation discussed in Section [4] Starting from the five-dimensional action S, we define 
a rescaled action S e which has the expansion 



Sf — Sn ~{~ € Sn 



(8.1) 



In Section [5] we found the most general solu tion of <5So = 0, and substituting that solution into S2, gave the four- 
dimensional action functional of Section |6 Cf . 



The basis of our approximation method is the smallncss of the bulk radius of curvature 1 /k n compared to the radius 
of curvature C c of the four-dimensional metric e x ^ a b- However, although this approximation is valid on all the branes, 



The action So for the solution is zero, assuming the brane tunings l|5.17|l 
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it inevitably breaks down as y — > ±00, far from the branes, as C c —> 0, as discussed in Section 4C It is worth noting 
that in the special case where all of the induced metrics on the branes are flat and there are no matter fields, the 
metric ansatz (with $ = const) is an exact solution to the five-dimensional Einstein equations, and this breakdown 
does not occur. 

One might expect contributions from the regime far from the branes to invalidate our four-dimensional effective 
description. However, we expect that the contribution to the action far from the brane will negligibly change the 
calculation, as in the region in which we expect large departures from the derived metric, the warp factor exponentially 
suppresses any contributions. 

It is possible for our two-lengthscale expansion to break down not only asymptotically, but also in between branes. 
A number of models (eg, [TU1 HH1 132] to cite but a few) discuss "bounce" behavior in the warp factor, where it 
decreases and increases again in between branes, as with a cosh 2 dependence. Typically, this behavior appears when 
the metric 7 is a curved FRW metric. It is a limitation of our method that this bounce is not evident in our solutions, 
as it explicitly requires coupling between the O(l) and 0(e 2 ) components (in particular, the four-dimensional Ricci 
scalar). Thus, this behavior is excluded by the underlying assumptions of our method, as near the turning point of 
these bounces, the separation of lengthscales has broken down. We note, however, that cosh 2 behavior is likely to be 
forbidden in the first or last (y —> ±00) regions by global hyperbolicity. It is also possible to produce sinh 2 behavior 
in the warp factor. In between branes, this can lead to topologically disconnected regions of spacetime as discussed 
in [28] . which we have excluded by assumption. In the first or last regions, correctly accounting for this behavior 
requires that the integration over the fifth dimension be truncated. However, the contributions to our effective action 
from integrating beyond these regions is again exponentially suppressed and negligible. In the regime in which the 
separation of lengthscales is valid, our solutions are in agreement with models displaying these types of behavior. 

For black holes, the solution given by our effective action is subject to the Gregory-Laflamme instability [33] and 
the final outcome is uncertain (see |34j and citations thereof). The five-dimensional stability of solutions for which the 
induced metric on the branes is not nearly flat (eg, neutron stars) is an interesting open question. We conjecture that 
all the solutions without horizons are stable and are reasonably described by our four-dimensional effective action. 

We may also consider the regime in which C c <C C, such as will occur a long way away from the branes. In this 
limit, the physical description would change from being that of decoupled fibres to that of decoupled four-dimensional 
hypersurfaces [one should solve the 0(e 2 ) contribution to the action first, and substitute that into the 0(1) contribution 
to the action]. This approach may yield a matched asymptotic expansion approach to obtaining a solution far from 
the branes. Our method may therefore be useful for investigating the regime between Minkowski space on a brane 
and a black hole on a brane. 

It is important to note that our method does not yield the leading order five-dimensional metric. This can be seen 
from the fact that our four dimensional action depends only on the fields "x evaluated on the branes, and the values of 
these fields between the branes are not determined. However, knowledge of the leading order five-dimensional metric is, 
rather surprisingly, not a prerequisite for correctly capturing the leading order four-dimensional dynamics. Most other 
methods rely on knowledge of the five-dimensional behavior of the metric to calculate the effective four-dimensional 
equations of motion, and our method is somewhat unique in this regard. 

Our method of computation correctly captures the leading order dynamics of the system. However, there will be 
higher order corrections, suppressed by powers of e 2 . In particular, the fields "x and ™$ can be expanded as 

n X = "X (0) + £ 2 "X (2) + 0(e 4 ), (8.2a) 

«$ = »$(<>) + e 2«$(2) + ( £ 4) (g 2b ) 

Throughout this paper, we have dealt only with the fields "x' '' and "<I >( - - > . The necessity for the higher order terms 
can be seen from the exact, five-dimensional equations of motion, which are derived in Appendix [A] For example, the 



exact Israel junction conditions are given by Eq. (A. 10). If we substitute the expansions ( 8.2[ ) into Eq. (A. 10), and 
use (5.7) [with n \ and ™$ replaced by "x^ and n $WJ together with the brane tuning conditions ( |5.17 ), we find that 
the higher order corrections n ^( 2 ' and n <3>( 2 ) are related to the matter stress energy tensors on the brane. Our results 
confirm the suggestion of Kanno and Soda that these higher order corrections do not affect the four-dimensional 
effective action to leading order [2"2"j . 



B. Models Which Violate the Brane Tension Tunings 



If a brane's tension is adjusted so as to violate the tuning condition (5.17), then it is possible to view the situation as 
having either detuned brane tensions or detuned bulk cosmological constants. For accounting purposes, it is simpler 
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to think of the bulk cosmological constants as being detuned. When this occurs, the exact equations of motion in 
the bulk (A. 6) imply that a nonzero Ricci curvature is induced to compensate for the detuning. Exact solutions have 
been calculated in highly symmetric cases, see for example Ref. [55]. In general, the exact nature of the perceived 
detuning is nontrivial, as the bulk cosmological constants on either side of the offending brane(s) can appear detuned 
by different amounts to compensate. 

If the deviation from the brane tuning conditions is small [Aa/a T = 0(e 2 )], then we can approximate the contri- 
bution to the four-dimensional effective action as 



jv-i . 

AS = - / d A x^/^h{a ri 



(8.3) 



where is the tuned value for the n th brane, given by (5.17). This approximation is of the same order as the other 



approximations we have made in our method. The net result is then an effective cosmological constant on each brane, 
given by 



Al 4 ) 



(8.4) 



which vanishes when the brane tensions are tuned. [Note: This differs from the result given in the literature for the 
RS-II model, see for example Ref. [5], but the difference is 0(e 4 )). 

If the detuning of a brane's tension from its tuned value should become too large [0(1) rather than 0(e 2 )], then 
the curvature induced by the four-dimensional effective cosmological constant can cause the radius of curvature on 
a slice of constant y close to the branes to violate the approximations used in our method, which implies that our 
four-dimensional effective action will not be a good description of a system in this regime. 



C. Multigravity 

Theories with more than one independent dynamical tensor field are called multigravity theories; see the general 
discussion in Damour and Kogan |29j . The models in this paper may exhibit two forms of multigravity, although we 
have ignored one of them entirely. 

The first form of multigravity is the possible existence of a second tensor field, given by the matrix B(x a ) in Eq. 



(5.10). We argued in Section 5C that this form of multigravity is likely forbidden. 

The second form of multigravity arises from the the fact that outside of the low energy regime, the models will 
contain Kaluza-Klein graviton modes. These modes will have masses that are formally of order but may be 

much lighter due to exponential suppression factors, and so be phenomenologically important (so-called "ultra-light 
modes" ) [25] [26] . Our method of analysis automatically excludes all massive fields (formally, we take e sufficiently 
small to overcome any large exponential factors), so we have neglected all graviton Kaluza-Klein modes. It is likely 
that some of these modes are in fact ultralight in our model, as in the analyses of Damour and Kogan [35J [5S] [53] • 



D. Conclusions 



The method we have developed is a useful tool for investigating the four-dimensional, low energy behavior of five- 
dimensional braneworld models. We have illustrated the method for a system of N branes in an uncompactified 
extra dimension. We intend to apply the method developed here to orbifolded models. We also plan to investigate 
experimental constraints on the model illustrated here, especially solar system constraints. Finally, we will investigate 
dark matter models obtained by placing matter fields on various branes. 
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Appendices 



Appendix A: Five-Dimensional Ricci Scalars and Exact Equations of Motion 



Here we present the dimensionally reduced Ricci scalar and the exact equations of motion for the action (3.24). We 
include the order at which terms appear in terms of our scaling parameter, e. 



Dimensional Reduction of the Ricci Scalar 



The constraint det-y = — 1 may be enforced either at the level of the equations of motion, or by using a Lagrange 
multiplier. 

If the constraint det-y = — 1 is being enforced at the level of the equations of motion, then it is simplest to compute 
the equations of motion using the metric (3.4), and then perform a conformal transformation on the quantities in the 



equations of motion. In this metric, the five-dimensional Ricci scalar is given by 

R (5) = £ 2 ( n R (4) _ 2V a Vq"$ \ _ n l ahn lab,yy + n ^ a hn„ 



1 9 3 

ra*2 V ' lab,y) > 



4"$ 



7 la 



1 ldb,y: 



(A.l) 



where covariant derivatives and the four-dimensional Ricci scalar are those associated with ™7 ah . 

For the constraint detj = — 1 to be enforced at the level of the action, a Lagrange multiplier term must be added 
to the action 



N 

AS = J2 d 5 x n n X(x a , y) (yf^ - l) , 



(A.2) 



where n \(x a ,y) are the Lagrange multiplier fields. Using the metric (4.5), the five-dimensional Ricci scalar is given 
by 



I / I n $ 
ncf,2 ( ~~A 7 7ac,y 7 ldb,y~ 5( X,y) ~ 4 X.yy + 4 — ^ X,y 



(A.3) 



where covariant derivatives and the four-dimensional Ricci scalar are those associated with "7 afc . To obtain this form, 
we use the following two formulae which may be derived from the fact that det(™7 ab ) = — 1: 



"I 11 " n iab,y = 0, 

n^ab n*ab fi£. n A/ c d n £j 

i !ab,yy / ibc.y I !da,y 



(A.4) 
(A.5) 



The complete action (with e scaling and Lagrange multipliers) is given by Eq. (5.1). 



2G 



2. Varying the Action 



We use ™7 ab to compute covariant derivatives, the four-dimensional Ricci scalar ( n R^) and the four-dimensional 
Einstein tensor ( n G^ b ). Indices will also be raised and lowered using this metric. 
Varying the action (3.24) with respect to ™$, we find the bulk equation of motion 

3 



:(V an x)(V "x)-3V a V a "x 



i(j)2 



Ky + ^2 ^ 7b c , a 7 7da, a ~ 2k 5 A « = °- ( A - 6 ) 



From combining the variations with respect to n 7 Qb and ™x (after eliminating the Lagrange multiplier by tracing 
over the ™7 afc equation of motion, or enforcing det 7 = — 1 on the equations of motion) , we obtain a traceless tensor 
equation of motion in the bulk 

\ n <f 2 e 2 e-^ U n G$ + n l ab n R (i) + 2(V a " X )(V b "x) - ^7 afc (V c "x)(V ' c n X ) - 4V a V b "x + n 7 ab V c V c "x) 



>$ e 2 e -"x (_4V a V b "$ +4(V (a "$)(V b) " X ) + " 7ah V c V c "$ - "7a 6 (V c "1>)(V c "x)) 



- n -y , -I !« n ^ , -2 n y n -y h + nA / n -y cdn ^y^ =0 

lab.yy 1 n ^ 'o-b.y A..y !ab,y 1 'ac,y ! 'db.y w 3 

and a scalar equation of motion in the bulk 



(A.7) 



+ 4 "7 afc "7 ah , ro + 3"x, ro + 3(«x,,) 2 - 3-^»x, w + 2"$ 2 . 2 A„ = 0. 



These variations also give rise to the boundary conditions on the branes 

1 „„ 



-■■vv _ n + 1 A, — 9 K 2 f 2 p~" x ( n T — "'V ^n~cdnrji 

$ lab,y n+ i $ lab,y ~ ZK 5 £ e I J ab 7ab 4 7 J c d J • 



and 



3 n X, u . 3" +1 x 



+ 



n+l(f> n O 



Here, the four-dimensional stress energy tensors on the branes { n T ab ) are defined by 



% S m [ n h ab + 6 n h ab , >] = n S m [ n h ab , n 4>]-^\ <£w n ^h n T ab 8 n h ah . 



(A. 



(A.9) 



(A.10) 



(All) 



Note that every factor of e 2 is accompanied by a factor of exp(— ™x). Also note that the 0(1) terms in these 
equations are exactly our equations of motion (5.3) to (5.7). 



